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Abstract 

We study the effects of a finite chemical potential on the occurrence of cavitation in a quark gluon 
plasma (QGP). We solve the evolution equations of second order viscous relativistic hydrodynamics 
using two equations of state. The first one was derived in lattice QCD and represents QGP at zero 
chemical potential and it was previously used in the study of cavitation. The second equation of 
state also comes from lattice QCD and is a recent parametrization of the QGP at finite chemical 
potential. We conclude that at finite chemical potential cavitation in the QGP occurs earlier than 
at zero chemical potential. We also consider transport coefficients from a holographic model of a 
non-conformal QGP at zero chemical potential. In this case cavitation does not occur. 
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I. INTRODUCTION 


The analysis of an extensive body of experimental information suggests that in the heavy 
ion collisions at high energies performed at the Relativistic Heavy Ion Collider (RHIC) [l|, l2| 


and at the Large Hadron Collider (LHC) j3(| a quark gluon plasma (QGP) is formed and 
it behaves like an almost perfect fluid with small dissipation. This highly energetic 


fluid is very well described by relativistic hydrodynamics js Sl. Some years ago it has been 
pointed out jo-lll that previously overlooked viscous effects could lead to processes such as 
cavitation in these systems. 

In general it is well known 12|] that if a non-relativistic liquid at a constant pressure 
reaches the critical temperature, the consequent ruptures of cohesion between molecules 
leads to boiling. On the other hand, cavitation occurs when at a constant temperature the 
critical pressure drops below the vapor pressure and vapor bubbles are formed. 

In an expanding quark gluon plasma, because of a non-vanishing viscosity, the pressure 
can become zero and even negative. In this situation cavitation would imply the creation 
of bubbles of a hadron gas phase in t he q uark gluon plasma phase. The evolution of these 
bubbles might have observable effects 13j. 

The study of cavitation hy the QGP was pioneered by the authors of Refs, js. 


further developed in Ref. [n|. Bulk viscosity has been a very active topic of research 
heavy ion physics [^. It is generally understood that for QCD at high temperatures, bulk 


101 and 


m 


viscosity effects are much sma 


limit is well approximated [m|. However, near the deconfinement region and in 


ler than shear ones, as for high temperatures the conformal 


ladronic 


15]. These 


systems, large values of the bulk contribution have been found in several works 
large values might trigger cavitation in the evolving plasma. 

In the present work we extend the formalism employed in [n| to the case of finite chemical 
potential. From phenomenological studies of particle production in heavy ion collisions with 
the help of thermal models we can conclude that the baryon chemical potential is of the 
order of /is ~ 1 — 10 MeV at LHC and pe — 20 — 100 MeV at RHIC. At lower energies, 
the chemical potential can reach higher values, fin — 500 MeV. Moreover, in the case of 
inhomogeneous freeze-out, even at high energies we can have domains of high chemical 


potential 


17l |. Recently the equation of state (EOS) of the QGP at hnite chemical potential 


was computed in lattice QCD simulations 18|, [l^ . In view of these results, we think that 
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it is interesting to investigate the effect of finite chemical potential on cavitation. In the 
search for different phenomenological parametrizations of the bulk viscosity coefficient and 
other transport parameters, based on different many-body micro scop ic dynamics, we have 


considered the holographic non-conformal model of QGP of Ref. 


20 | 


II. SECOND ORDER RELATIVISTIC VISCOUS HYDRODYNAMICS 


A. The main equations 


Throughout this study we use natural units h = c = ks = I- The energy momentum 


tensor of viscous relativistic hydrodynamics 


[n, 13, 20 


22 j is: 


= eu^^u'^ - 


( 1 ) 


where the e and p are the energy density and the pressure of the fluid, = ( 7 , yff) is 
the 4-velocity of the fluid (7 is the Lorentz factor 7 = (1 — with normalization 

= —1. The metric is given by g^i, = diag(—,-f,-f,-|-). fl^'^ is the viscous tensor 
satisfying = 0 and + u^u'^ is the projector orthogonal to u^. Notice that 

we do not include dissipative currents associated with finite chemical potential and heat 
flow. In Bjorken flo w l24j| . symmetry constraints eliminate these contributions from the 
equations of motion 2^. The heat flux satisfies the orthogonality property = 0 in 
the Landau-Lifshitz frame [^. The two relevant terms in the coupling of the dissipative 
equations involve the terms V<^q'i/> and In the Bjorken flow the temporal component 

of the heat flux is zero {q'^ = Oj^which leads us to the conclusion that the previous heat flux 


contributions are exactly zero 


25| 


Conservation of the energy and momentum is enforced by = 0, where is the 

covariant derivative. The projection of this equation on the direction parallel to the fluid 
velocity reads: = 0. Similarly, the projection on the direction perpendicular to the 


fluid velocity 


equations 


and 


11 


vie 


20 


s: = 0. These two projections lead to the following evolution 


22 |: 


De -\- [e p')9 V_L {^Uu) — 0 

{e + p)Du^ + v^p + A^V^n“^ = 0 


( 2 ) 


(3) 
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In the above equations we have D = u^Vfj,, 9 = = A^"Vq and the symmetrization 

of the indices is represented by + Ai,^)/2. 


The viscous tensor is given by [11 


,Q,Q: 


W'' = 


(4) 


where is the traceless symmetric part tt^ = 0 and 11 = We obtain the 

components of the viscous tensor® solving the following equations [ul, l20|, l22| which are 
second order Israel-Stewart like |23|: 


r, 


n 




n = -C 0 - 4T)n 


(5) 


( 6 ) 


with the definitions A^^“^ = (A^^A*^^ -|- A^^A^"") /2 — A^^A“^/3 and 


The shear viscosity is rj, and the seconc 


holographic calculations of TV = 4 SYM 26 


order coefficients are and Ai, given by the 


27|: 


r] 1 

s dvr 


Yi “ 


[2 - /n(2)] 


27rT 


Ai = 


r] 


27rT ^2^- 

and by the following parametrization of C/s [ll|, |21| : 


(7) 

( 8 ) 
(9) 


C [t^-t 

s AT 




for T>Tr 


( 10 ) 


whith a = 0.901, b = 0.061 and AT = Tc/14.5 . We also consider = r^. We obtain the 
entropy density from the relation 


s = 


dp 


( 11 ) 


For this parametrization we will choose Tc = 190 ~ 200 Mel/. Setting the second-order 
coefficients equal to zero (r^ = 0, = 0 and Ai = 0) we recover the Navier-Stokes approach 

in (jSj) and ([6]): = —2pa^'' and 11 = —(,9. 
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B. Bjorken flow 


Now we consider solutions of the hydrodynamical equations which do not depend on 
the transverse spatial coordinates x and y and are boost invariant in the 2 ;—direction as 


in 


11 


20 


2l|, which are the Bjorken symmetries. To obtain the evolution equations in 

rom {t,z) to (r, 0 , 


this boost invariant description we perform the changes of variables 
where r is the proper time and ^ is the spacetime rapidity, given by 2 ^: r = y/t^ — and 
^ {z/1). In general, the Milne coordinates are (r, r, 0, and the line element is 

ds^ = —+ dr"^ + so 5 ^^,^ = diag{—l, 1,1,r^). Due to the Bjorken symmetry 

the four-velocity in Milne coordinates turns out to be = (1, 0, 0, 0). 

In the (r, coordinates we have 9 = 1/r, D = Vr, tL = tL = 1/r, FT = r and 


= —T^j^ux = Tjjj, . In particular, = r . 
The energy-momentum tensor ([1]) becomes: 

/e000\ /o 0 


= 


0 p 0 0 
0 0 p 0 

y 0 0 0 p y 


+ 


0 n +1 


0 
0 
n + 


0 

0 

0 


\ 


( 12 ) 


0 0 

\o 0 0 n-<i>y 

where the contributions to the pressure from the bulk and shear stresses are respectively 
given by If (the trace of 11"^) and $ (the traceless part of 11"^), given by = <h/2 

and 7r| = — <h. 

From (1T2|1 we identify the transverse pressure: 

= p + n 

and the longitudinal pressure: 

Fg = p -I- n — $ 

In tl^ (f 0 description, the equations ([2]), ([3]) and 
tions 


$ 

'2 


Q: 


(13) 


(14) 

give the following evolution equa¬ 


te 


£ -I- p -t- n — <F 


,a<F 


and 




an 


4tS 


3r 


2 p2 


^1711 _ 


or T 


(15) 

(16) 

( 17 ) 


Cavitation will take place when the longitudinal pressure becomes negative. 
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III. EQUATION OF STATE 


In this work we study the occurrence of cavitation considering two equations of state for 
the quark gluon plasma. 


A. Model 1 

The following EOS was previously employed in [n| and we use it to perform comparison 
with other models. This equation of state describes the quark gluon plasma phase and the 
crossover to a hadron gas, being a parametrization of a lattice calculation at zero chemical 


potential with Tc = 190 MeV. The trace anomaly is given by 

/ \ 


- 3pi 
2^4 


V 


1 + exp 


T-ci 

C2 


1 2 


/ 


^2 ^4 


(18) 


d? = 0.24 GeV^, = 0.0054 GeV^, Ci = 0.2073 GeV and C 2 = 0.0172 GeV. The pressure 

is: 

r / 

1 





1 - 


1 + exp 


T -Cl 
C2 


/ 


d2 d^ 

J^'3 5 


(19) 


From firsj) and (IT^ we obtain the energy density: 

/ 

ei=T^ 


1 - 


[i+(kf)]) 


^2 di 
rp 2 ' 2^4 





1 + exp 


T'-ci 

C2 


^2 di 

2"'3 J''5 


( 20 ) 


B. Model 2 


In this work we study finite chemical potential effects in cavitation. To do so, we use the 
recent parametrization of a lattice simulation of SU{3) QCD matter at finite temperature 
and chemical potential, with quarks {u, d and s with equal masses) and gluons 18l. Il9|. 

We have previously used this equation of state in the study of the expansion of the early 


Universe in 


29| . where we have developed the expressions in detail. We start by recalling 
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the trace anomaly at finite chemical potential 18|, Il9l |: 

P2{T,fi) 


g2(T,/i) - 3p2(T,/i) ^ 


2^4 


I ^ 

+ J^X2 


^ £2(T,0) - 3p2(T,0) dX2 

T4 ^ 2T dT 

where the chemical potential contribntion is given by the function [l 8 |: 


X2{T) = • /3 • 


In the zero chemical potential limit we have 


18 


tanh{fi • r + /g) + 1 


( 21 ) 


( 22 ) 


19|: 


£ 2 (^, 0 ) - 3p2(T, 0) 

2^4 




hn + 


/o • 


tanh{fi ■ r + / 2 ) + 1 


1 + Qi-T + g2-T‘^ 


(23) 


In the last three expressions we have introduced the variable r = T/200 MeV, where 200 
MeV is the critical temperature. 


We consic 
(i) From 


er two sets for the dimensionless parameters: 


lOj: ho = 0.1396, h = -0.1800, ha = 0.0350, /o = 1.05, /i = 6.39, /a = -4.72, 


gi = —0.92 and 22 = 0.57 . For the chemical potential parametrization we use from 
hg = -0.5022^4 = 0.5950, /g = 0.1359, U = 6-3290 and /s = -4.8303 . 


18|: 


(zi) From ho = 0.1396, hi = -0.1800, /^= 0.0350, /o = 2.76, /i = 6.79, /a = -5.29, 


2 i = —0.47 and 22 = 1-04 . This set is used in 
We also consider this set in our section IV. 
The pressure is calculated from fl2ll) : 

*r 


20 [ | for the non-conformal holographic plasma. 


Pa(T,p) = T' 


dT 


, —hi/T 


a 


r 


hn + 


/o • tanh{fi ■ r' + /a) + 1 


1 + gi-r' + g 2 -T' 


+ (24) 


Inserting (124)) into (| 2 T]) we find the following expression for the energy density 29| : 

/o • [fanh(/i ■ r + /a) + 1 


e 2 (r,p) ■ 


^ ^-hx/r'-hi/r' 


hn + 



r 


hn + 


1 + 21 • r + 22 • 


/o • tanh{fi • r' + /a) + 1 


2 dT 

(25) 


1 + 2 i • 't' + 22 • 

In Fig. [T] we present the trace anomaly given by {e — 3p )for Models 1 and 2 at zero 
chemical potential and for Model 2 at finite chemical potential. 
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(a) 


FIG. 1: Trace anomaly for several chemical potentials. 


C. Numerical Results 

We start our numerical study following the calculations developed in [h|, with no bulk 
effects and considering that the evolution starts at tq = 0.5/m, with d>(ro) = 0. We want 
to compare different equations of state and determine the effects of certain features of these 
EOS (such as, for example, finite chemical potential) on the evolution of the fluid. To do 
this we would like to evolve them starting from the same initial state. However this is 
not possible because if they represent systems at the same initial temperature, then, these 
systems will have different initial energy densities and vice-versa. In relativistic heavy ion 
collisions, before thermal equilibrium formation, energy is released from the projectiles in 
a certain volume in the central rapidity region. It is perhaps more realistic to say that the 
quark-gluon system has first a defined energy density determined from the conditions of 
the collisions and then it reaches thermal equilibrium, forming a QGP with the properties 
determined by the equation of state. Since it is very difficult to know what comes first, 
energy density or temperature, we will test the two possibilities: same initial temperature, 
Tinitiai = 305 MeH and same initial energy density, Einuiai = IQGeV/fm^. We will find out 


that our conclusions regarding cavitation do not depend on this choice and are the same for 
different initial conditions. 








We solve the coupled equations (IT^ and flTdl) to obtain the temperature for each model 
considered. The temperature evolution at zero chemical potential is showed in Fig. [2l The 




X (fm) 

(b) 

FIG. 2: Temperature evolution with no bulk effects for Model 1 and Model 2 (i) from tq = 0.5 fm 
to Tfinai = 10/m at zero chemical potential, a) Same initial temperature, b) Same initial energy 
density. 

effects of the chemical potential are presented in Fig. 12] for Model 2. From the hgure we 
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clearly see that systems with larger chemical potential cool faster. 



X (fm) 


(a) 



T (fm) 


(b) 


FIG. 3: Temperature evolution with no bulk effects for Model 2 (i) from tq = 0.5 fm to Tfi^ai = 
10 fm at finite chemical potential, a) Same initial temperature, b) Same initial energy density. 


Cavitation is an effect mostly associated with peaks of the bulk viscosity close to phase 


transition 




We now consider bulk effects via eq. (HU with the relaxation scale 
approximated by = r^. As previously mentioned, cavitation is generated when the 
longitudinal pressure, P^, is negative and it starts when P^{Tcav) = 0 in flT^ . where Tcav is 
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the time at which the longitudinal pressure vanishes. In this subsection we will make use 
of d?]) to flTT]) to solve numerically the three coupled evolution equations ([15]), flT^ and flT7|) 
and determine the time evolution of the longitudinal pressure fll4p . 

We start with zero chemical potential in Fig. 0] for the two models with the initial 
conditions Tq = 0.5/m, <l>(ro) = 0, n(ro) = 0 for Ti^uiai = 305 MeV and also for Einuiai = 
16 GeV/ fm^. 


TABLE 1: Time when cavitation starts (in Fig. 4(a) 


Tcav [fm] 

(fixed Tinitiai = 305 MeV) 

Model 1 

2.3 

Model 2 (i) 

2.9 

TABLE 11: Time when cavitation starts (in Fig. |4(b)|) 

Tcav [fm] 

(fixed Einitiai = lOGeV/fmf) 

Model 1 

2.22 

Model 2 (i) 

4.06 


In Fig. [5] we show results at finite chemical potentials for the Model 2. The cavitation 
starts earlier for increasing the chemical potential. 


TABLE III: Time when cavitation starts (in Fig. 5(a) ) 


Tcav [fm] 

(fixed Tinitiai = 305 MeF) 

Model 2 (i) (p = bOMeV) 

2.91 

Model 2 (i) (p = 250 MeV) 

2.84 

Model 2 (i) (p = 500 MeV) 

2.63 


IV. CAVITATION IN A NON-CONFORMAL HOLOGRAPHIC PLASMA 

One interesting question is whether it is possible to hnd cavitation within different ap¬ 
proaches for the underlying many-body microscopic dynamics. From the strongly coupled 
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(b) 


FIG. 4: Time evolution of the longitudinal pressure (I14p from tq = 0.5 fm and zero chemical 
potential. See Tables |T] and |TT1 a) Same initial temperature, b) Same initial energy density. 


nature of the QGP, in particular near the deconfinement transition temperature, the micro¬ 
scopic physics should be inherently non-perturbative. In this context, holographic methods 
may shed light upon the nature of transport coefficients for such strongly coupled plasmas, 
since these coefficients are difficult to obtain from lattice methods. 

We follow the parametrization and second-order theory from 20|. Considering Bjorken 
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(a) 



(b) 

FIG. 5: Time evolution of the longitudinal pressure (jl4|) from tq = 0.5 fm and finite chemical 
potential for the Model 2. See Tables IIIII and IIVI for values of the chemical potential, a) Same 
initial temperature, b) Same initial energy density. 

symmetry and flat spacetime, our analysis reduces to solving the following dissipative equa¬ 
tions: 

^ ( 26 ) 

\ / V 'jC 
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TABLE IV: Time when cavitation starts (in Fig. |5(b)| ) 


Tcav [fm] 

(fixed Einitiai = \0GeV/frr?) 

Model 2 (i) (p = 50 MeV) 

4.04 

Model 2 (i) (p = 250 MeF) 

3.8 

Model 2 (i) (p = 500 MeF) 

3.15 


and 


{DU + n^) + n = -ce + Din ( ^ 




The employed holographic method is a bottom-up solution 


(27) 


30| that parametrizes some of the 


thermodynamical properties of lattice (2-1-1)-flavor QCD near the crossover phase transition, 
described by Model 2 (ii) aX zero chemical potential with the parameters [l8|, |20|. The trace 
anomaly with these parameters is shown in Fig. [H 


It is known that ho 
value of r]/s = l/dvr 


ographic calculations of transport coefficients, such as the holographic 


26| . are generally consistent with experimental investigations of the 


evolution of QGP. The calculations employing usual kinetic theory result in larger values of 
7]/s, which corroborates that holographic calculations are at least consistent with the QGP 
evolution near the deconhnement region 31|. 

The advantage of this approach is that the non-trivial temperature dependence of 
the transport coefficients were calculated for a phenomenological bottom-up model and 
parametrized in 20|. In the next subsection we summarize the relevant results. 


A. Transport coefficients and Bjorken flow 


The coefficient is still l/dvr, as is expected for several holographic models 
The shear relaxation time is parametrized as follows {T^ = 143.8 MeV): 


26 


32| 


v 

^nV 


T \ 


X = — \ = 


T2 J'h J ^ ^bh{ch-x) ^d^{e^-x) ^fh{gh-^) ’ 

where to are fit parameters listed in Tab. |Vl 
The parametrization of bulk viscosity is: 

T 


X = 


'rh 

^ C 


^bulk _j_ u'utft: 

btikT + {cLikf 


jh 

^bulk 


( 28 ) 


(29) 
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where to e^if. are fit parameters in Tab. |Vl 

In order to parametrize the bulk relaxation time, we adopt a causality bound, relating it 
to the shear relaxation time 20 | : 


4t{x = 






+ (cj) 


rjh 

=+x 

2 X ' 


(30) 


with the corresponding £t parameters to being given in Tab. IVl 


TABLE V: Parameters: r^rj/T^ for Eq. (I28]l . (/s for Eq. (f2^ and r^T for Eq. pOl) . 



bh 

V 


Jh 

^T) 


fh 

•>v 

4 

0.2664 

2.029 

0.7413 

0.1717 

-10.76 

9.763 

1.074 

„h 

%ulk 

bh 

^bulk 

Ji 

'-bulk 

'^bulk 

^bulk 



0.01162 

1.104 

0.2387 

-0.1081 

4.870 





„h 

d'l 




0.05298 

1.131 

0.3958 

-0.05060 





The other coefficients are determined in a phenomenological approach. However there is 
still a non-trivial behavior near the transition temperature. These are estimates considering 


fluid-gravity ca 
calculations 27 


culations for conformal fluid, as well as other strongly coupled holographic 


33 


34j |. We consider then the following 20| : 


^ sT ’ 


T* — — 

Tr — 'Jpl 


6 = Ai 
6 = ‘^VTn 


- c. 


- C, 


- C, 


(31) 


In Fig. 6 (a) we show the two bulk expressions given by (fTUD and 


parametrizations from 


[n|, y, 


using the 


21| . Notice that near the deconfinement temperature the 
bulk viscosity per entropy density of the holographic calculation given by fl29]) is significantly 
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smaller than Both come from distinct microscopic approaches and the difference be¬ 

tween them reflects our lack of knowledge on the subject. We hope that future theoretical 
and experimental work can shed light upon bulk viscosity properties of strongly coupled 


plasmas. Similarly, in Fig. 6(b) we present the two relaxation times as a function of the 
and 


temperature, (18J) and fl28]) . Eq.(|H]) is the result of a conformal calculation and Eq. 
includes non-conformal contributions calculate through holography. 

It is straightforward to solve Eqs. fl26|) and fl27|) in the Bjorken flow. We follow the 
conventions adopted in the previous section, so the evolution equations are flT5|) and the 
following: 


"ar 3 t 


-^1 . t .2 


4r^ 

3r 2r]^ 


+ 

3 C 


and 


an 


a c^rr ^ XT C. T-r 3 ^1 9 ^2 9 

^-§7 + ^n 7 =- n +-hi’''+ htf + 


4 


A 

'WT 


(32) 


dT 

(33) 


2r/2^ ' ' (C/s) 

In particular, equations fl3^ and fl3^ are the same as ffT6|) and flTTI) but with extra transport 
coefficients. One interesting feature of such second order theory is that both bulk and shear 
evolution equations are nonlinear and coupled to each other. The entropy density in fl2^ . 
fl32ll and fl33ll is calculated via flTTll using the equation of state of the Model 2. 

We solve numerically Eqs. flT^ . fl3^ and fl3^ in the Bjorken flow and with the discussed 
holographic transport coefficients to search possible indications of cavitation. We consider 
the initial conditions Tq = 0.5 fm, <F(ro) = 0, n(ro) = 0 for both initial approaches Tinuiai = 
305 MeV and also Sininai = IQ GeV/ fm^. 

As can be seen in Fig. 0 cavitation does not occur. It is often the case with holographic 
transport coefficients that dissipative effects seem to be smaller than others obtained from 
different methods, such as kinetic theory calculations. Nonetheless, the holographic values 
for small ratio of shear viscosity per entropy density are consistent so far with experimen¬ 
tal results, and thus we suggest that, if conhrmed by experimental data, the absence of 
cavitation may be a consequence of the strongly coupled nature of QGP. 

In order to compare the Models 1 and 2 with the non-conformal holographic Model we 
show again, in Fig. |H] and in Fig. [9] the same plots of the temperature and longitudinal 
pressure evolution previously presented in Fig. El Fig. |4]and Fig. [71 


16 










T (MeV) 

(a) 



(b) 


^IG. 6: fa) Bulk per entropy density, (b) Relaxation scale. The parametrizations come from 

111, l2o|,l2li. 


V. CONCLUSIONS 


In this work, we considered several equations of state at different chemical potentials in 
order to investigate the occurrence of cavitation in QGP. This was done through a numerical 
integration of the Bjorken hydrodynamical equations and determination of the pressure as 
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T (fm) 

(a) 


3 


(GeV/fm") 


Non-Conformal Holographic QGP 


•T , = 305 MeV 

initial 

e . = 16GeV/fm’ 

initial 


0.5 2.0 3.5 5.0 6.5 8.0 9.5 

T (fm) 


(b) 


FIG. 7: (a) Temperature evolution for the non-conformal holographic QGP. (b) Time evolution of 
the longitudinal pressure ()14l) for the non-conformal holographic QGP at zero chemical potential. 


a function of proper time. 

We used a phenomenological parametrization of the lattice QCD equation of state with 
two sets of transport coefficients. For the bulk viscosity we used a parametrization of results 
obtained in a Monte-Carlo simulation of pure gluon dynamics SU(3) UJ, l3^. For both 
lattice QCD and mean held QCD equations of state, cavitation occurs at smaller time scales 
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T (fm) 

(a) 



T (fm) 

(b) 

FIG. 8: Temperature evolution with (a) same initial temperature and (b) same initial energy 
density. 

as the chemical potential increases. 

As is often the case with transport coefficients calculated with lattice methods near 
the phase transition, the height of the bulk per entropy density peak is not generally well 
established, most calculations rely on pure glue dynamics and consequently on the existence 
of hrst order phase transitions. We considered a phenomenological holographic calculation 
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(b) 


FIG. 9; Longitudinal pressure evolution with (a) same initial temperature and (b) same initial 
energy density. 


of transport coefficients of a non-confornial QGP 20|. In our approach, we conclude that 
cavitation does not occur in the holographic prescription, specially due to the small bulk 

). Even with the addition of more transport 


viscosity as shown in Fig. 6(a) and by Eq. 
coefficients (which contain peaks near the crossover phase transition) this approach does not 
support cavitation. In the future it may be possible to compare these distinct cavitation 
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results in the scenario where lattice calculations of bulk effects can consistently consider the 
crossover phase transition. We hope that once the magnitude of bulk effects in the QGP is 
more precisely calculated, we can be more conclusive on whether cavitation occurs. 

It should be interesting to generalize our calculations to more complex flows and check 
whether cavitation can be obtained in different regimes. Also, a non trivial check of the 
observable effects of cavitation in QGP would be the study of stability properties of hadron 
gas bubbles, as well as how their collapse influences the overall evolution of the plasma. We 
leave such investigations for future work. 
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